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Abstract. We continue the study, from a semiclassical viewpoint, of Calvo and Picon's 
operators, as manipulating photon states in quantum communication. In a previous 
paper, we defined a one-dimensional model operator and studied it analytically before 
moving on to Calvo and Picon's full two-dimensional operators. In the present paper, 
we show how the one-dimensional operator may also be useful as an experimental model, 
since it allows manipulations of two-dimensional Laguerre-Gaussian modes; the intensity 
distributions (in physical space) of the Laguerre-Gaussian modes then approximately 
flow along the elliptic curves studied earlier. Since the Wigner transform is fundamental 
in the study of Laguerre-Gaussian modes, we give a slightly expanded and improved 
treatment of the semiclassical Wigner transform, which was only briefly mentioned in 
the previous paper. 



I. Introduction 

Gabriel F. Calvo and Antonio Picon [3j introduced a class of operators which allows 
arbitrary manipulations of the three lowest-order Hermite-Gaussian (HG) modes for the 
purposes of quantum communication. In our previous paper [6], we (1) classified the 
self-adjoint extensions of the generators and (2) showed how to construct semiclassical 
approximations of the associated unitary operators; for the latter, we first considered a 
simple model operator to illustrate the semiclassical methods in a simplified setting before 
then moving on to the full operators of Calvo and Picon. 

In this paper, we show that the simplified operator also serves as a model operator 
from an experimental viewpoint, when considered as manipulating Laguerre-Gaussian 
(LG) modes of orders (0,0), (1,0), (0,1), and (1,1). 

Calvo and Picon introduced the following eight generators acting on Hermite-Gaussian 
modes: 
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defined in terms of tfie creation and anniliilation operators dl. = (^^ ~ ^) and dx = 
^ (a: + ^) , respectively, and similarly for the variable. 

These generators, within the subspace generated by the lowest three Hermite-Gaussian 
modes Hf = {|0, 0), |1, 0), |0, 1)}, obey the SU{3) algebra 

(a, 6, c = 1, 2, ... 8), where the only nonvanishing (up to permutations) structure constants 
fabc are given by 

/l23 = 1, /l47 = /l65 = /246 = /257 = /345 = /376 = 1/2, and /458 = /qts = 

We note that the triad of generators 

Ti^ {71,7-2,7-3} 

gives a SU{2) group that conserves the mode order. The remaining two SU{2) groups 
are formed by the triads 

r2 = {%,%,{f3 + V3fs)/2} 



and Ts = {%, fj, (-7-3 + V3fs)/2}. 



Unitary operators Un generated by the first triad give rise to superpositions between the 
two modes |1, 0) and |0, 1), leaving invariant the fundamental mode |0, 0). Unitarities 
and Ur^, generated by the second and third triads, produce superpositions between the 
two modes |0, 0) and |1, 0) (leaving invariant |0, 1)), or the modes |0, 0) and |0, 1) (leaving 
invariant |1,0)), respectively. 

Here we instead consider the one-dimensional operator 



T = d^ + d 



d^d^d 



a^dd. 



acting on the one-dimensional Hermite functions This is a simplified version of the 
operator above. We introduced this operator in |6] since it exhibits the same essential 
behavior as 7^, while leaving out the unessential additional variables. Here we will show 
that this operator may also be useful in experimental work, as an easier preliminary step 
before moving to the more complicated operator 74. 

One may check that Th^ = hi and Thi = Hq. More generally. 



Thr 



where 



(3n = {l~n)Vn + l. 

So T, when acting on the domain of (finite) linear combinations of HG modes, behaves 
precisely as the infinite Jacobi matrix 
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Using a theorem from Berezanskii's book ([2], p. 507; see also [6|), one can show that 
the deficiency index of this operator is (1,1), and hence the space of boundary values is 
of dimension 2 = 1 + 1. Moreover, using results of Allahverdiev [1], one can explicitly 
classify all self-adjoint extensions in terms of certain boundary values at infinity. This 
was achieved for the full operators of Calvo and Picon in [6], using a slight modification 
of Allahverdiev 's methods, but for the simplified operator T the methods of Allahverdiev 
may be applied without any modification. 

We will show that this model operator T allows simple manipulations of LG modes. 
But first, since the Wigner transform plays a basic role in the theory of LG modes, we 
recall some relevant facts in the next section. 



II. The Semiclassical Wigner Transform 
The standard n-dimensional semiclassical Wigner transform is defined as 

(1) W{f, g) (x, = (27r/^)-"/2 j e~^^^/'f{x + ^p)g{x - ^p) dp, 

where / and g are square-integrable functions of n variables. Then, among other useful 
properties (one may consult Folland's book ^), we have the norm-preserving property: 

\\W{f,g)\\L2(R2u) = ||/||j^2(Kn)||5(||i2(JRn). 

Also, the Wigner transform has an important relationship with Weyl quantization. We 
recall that the semiclassical Weyl quantization of a symbol a is the operator Op^(cr) 
defined by 

OpY{cT)f{x) = i^^h)-" 11 e''^^-y^^/'^a{{x + y)/2,0f{y)dydt 

Then 

(Opr(a)/|(7) = {2nh)--/' J j' a{x, OW{f, g){x, ^dx d^. 

As discussed in the previous paper p], we have an approximate formula for the evolution 
of the Wigner transform. Let X be a self-adjoint semiclassical pseudo differential operator 
with Hamilton flow nt generated by the (possibly /i-dependent) Weyl symbol, and let 
Ut = e^**'^/'* be the unitary propagator. (Note: In the previous paper, Ut sometimes 
denoted the semiclassical approximate unitary propagator.) We then have 

iy(/,^7)(«:_i(x,0) = W^(?7t/,?7i(7)(a;,O + 0(/i')||/||L2(M")||^7||LW. 

This essentially follows from Egorov's theorem with an 0{h'^) error (see [6] and references 
therein) . 

In the previous paper, we briefly discussed two-dimensional Wigner transforms. In the 
present paper, in the next section, we will consider the case n = 1, since LG modes are 
themselves one-dimensional Wigner transforms of Hermite functions [S]. 
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III. The Semiclassical Laguerre- Gaussian Modes 

In the context of LG modes, it is more convenient to use an alternative definition of 
the Wigner transform. We define the (one-dimensional) extended Wigner transform, of a 
function F of two variables, as 



iy(F)(x,0 = (27r/i) 



(The sign of the phase is intentionally different from ([T]).) This is clearly a unitary 
operator; moreover, as shown in [5], the (m,n)*'* LG mode is precisely given by Whmn, 
where hmn is the {m,ny^ HG mode. For reference, we will now state the main facts in 
the semiclassical setting. 

The zeroth order HG mode is the Gaussian function 

/ioo(a;,y) = (7r/i)-i/2g-(.2+,2)/(2h)_ 

All other HG modes, hmn, may be recovered by applying the creation operators 
a\ = {2h)-'/^ (^x - and dl = {2h)-'/^ (^y - 

That is, 

The corresponding annihilation operators are 



ai 



{2h)-'/^ (^x + and d2 = {2hY^'^ + 



Alternatively, to recover the LG modes, one applies to /iqo the creation operators 

A^, = —;=(a\ + iai) and AI_ = —;=(a\ — ial). 
^ V2 V2 

The corresponding LG mode annihilation operators are 

y4+ = —={ai — and A_ = —={ai + ^2). 

V 2 V 2 

As discussed in [5] (but also as one may easily check), the extended Wigner transform 
intertwines the two classes of creation and annihilation operators: 

(2) ^ . . r . . 

A+W = Wdl, A_W = Wd2. 

Since WhoQ = /iqo, this shows that the LG modes are indeed the extended Wigner 
transforms of the HG modes. Moreover, one may deduce the following familiar analytic 
expressions for the LG modes. Suppose x, ?/ G M, and let ^ = x + iy. Then, with hjk 
denoting the (j, kY^ HG mode, the [j, ky^ LG mode is 

~ J(7r/i)-i/2(A;!/j!)V2(_i)^'(^/v/^)i-fee--/(2^)Lf '=(2z//i) if j > k, and 



where 



L^(x) = ^_^-^(e-^a;"+") 
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are the Laguerre polynomials. Of particular importance in this paper are L'^ix) = Lq[x) = 
1 and Ll{x) = 1 — x. 



IV. Manipulation of Semiclassical LG Modes 
We first consider the action of the tensor product T®T on tensor products of functions: 



{r^r){f^g) = T{f)®T{g). 

Then we simply pull back this action to operate on extended Wigner transforms of func- 
tions. We write this as follows: 



irW)if0g) ■.= W{T{f)®T{g)). 

Since the extended Wigner transform has the intertwining property ([21), we see that the 
pullback operator T* is precisely the differential operator 

{a\ + i+ - a\a\a+ - a\a+a+) o (it + i_ _ aIaIa^ - aIa^a^). 

However, this point of view only seems to complicate matters. 

We are especially interested in the case when f ®'g = hmn, the {m^nY^ HG mode; then 
the pullback T* acts on LG modes. In particular, T acts on the four "binary" HG 
modes as follows: 

(T(g) T)/ioo = /ill, (T® T)/iii = /loo, 

(t ® f )hio = /loi, (t ® t)/loi = ho. 

Hence 

{f*W)hoo = Wh,, {f*W)h, = Whoo, 

(t*iy)/iio = Whoi, {f'W)hoi = Who. 
These lowest-order LG modes may be written in polar coordinates as 

Whoo = (7r/i)-^/2g-^'/(2h)^ ^f^^^ ^ -{nh)-'/\l - {rlVhf)e-'-"/^^^\ 

Who = (7r/i)-i/2(r/v^)e*^e-'^'/(2/.)^ ^^^^ ^ {^h)-^'\r / Vh)e~'' e~'-" /^^""l 

So we see that the pullback T* interchanges the LG modes Whoo and Whi and simply 
causes a phase reversal in the LG modes of order one. 

For convenience, in order to have a semiclassical differential operator, we define 

T = -2-^i^h'^f, 

and we note thalEI 



% = + (hDf) - 2hx - lihD, 



with the notation D = 

I ax 



^In [6] we mistakenly studied the operator \x{x^ + (hD)^) ~ — \ih^D with Weyl symbol \x{x'^ + 
— which is not quite the same as T. However, it was merely given as an example, to illustrate a 
general method, so in that paper one may simply redefine the operator in the first place. 
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Figure 1. For convenience we let /i = 1 and T = t/\/2, so that 
W^^WhoQ = 7r-i/2e-(^'+2/')/2[cos(2T) - ysm{2T) + (x^ + y^) sm\T)]. The 
figure shows the absolute value for T = /c7r/8, k = 0, . . . ,S. 



Just as we defined the puUback of T, we may similarly define the puUback of the unitary 
propagator Ut = e"**'^/'^: 

{U:W)if0g) = WiUtf®U;^). 
One may check that Ut acts as follows: 

Utho = cos{{h/2y/H)ho + ism{{h/2y/H)hi. 

Hence we have 

U*Whoo = cos'^{{h/2y/H)Whoo + sm'^{{h/2y/H)Whu 

+ I cos{{h/2f'h) sm{{h/2f/H)[Whro - W^V], 

where, in polar coordinates, 

Whio - Whoi = {-nhy^/^r / Vh)e-''''^^''\2isme). 
This transformation is pictured in Figure [H 
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As mentioned in Section [TTl we have an approximate formula for the evolution of the 
Wigner transform in terms of the Hamilton flow. At the end of this paper, we will 
adapt the result to the case of the extended Wigner transform, for the study of LG mode 
intensities; but first we will explicitly compute the Hamilton flow. 



The Weyl symbol of T is 



p{x,^; h) = + ^'^) — 2hx. 



Later we will need a re-scaled version of p, so we consider the slightly more general symbol 

Pr{x, ^; h) = ]^x{x^ + - ^r^/ix. 
Hamilton's equations for this symbol are 

X = x^ 



(3) 



with the conserved quantity C = pr{x,C,] h). 
First suppose that C ^ 0. Then we have 

a;{t) = -c(pit + to)--r^h] ' and - "^^^ + ^^^^ 



where to is either an arbitrary real constant or an arbitrary real constant plus the 
purely imaginary half-period of p. Here p is the Weierstrass p-function associated to the 
invariants 

= Lr^h^ and ^3 = -C^ - —r'^. 
^ 12 4 216 

For further details, one may consult our previous paper [6]. 

When r = and C 7^ 0, ^(t) is always strictly decreasing, which follows simply from 
Hamilton's equations Q. However, when r > we have a more complicated behavior, as 
shown in Figure [21 There is a pocket of radius r^/h = :^u7o, where wq is in practice the 
radius of the laser beam's waist. 

For C = 0, depending on the initial conditions, we have one of the following four cases: 

Vt 




= v^coth (^(t + to)) , 
'x{t) = Vt 

^(t) = v^tanh (^^(t + to)) , 

x{t) = ±\^r'^h sech {\/rVi{t + to)) 
^(t) = -Vr^tanh (yrVi{t + to)) , 
or 

|(x(t),^(t)) = (0, ±Vr^h) Vt (hyperbolic stationary points). 

Of course, we also have the elliptic stationary points (x, ^) = (±-y/ 0), corresponding 
to C = T{r^h/3f/^. 



8 



MICHAEL VANVALKENBURGH 




Since the LG modes are extended Wigner transforms of HG modes, we now see that 
the intensity of LG modes evolves according to the Hamilton flow we have just computed. 
As discussed in Section [Tll we have that the evolution of the Wigner transform may 
be approximately described by the Hamilton flow Hf However, the extended Wigner 
transform acts slightly differently with respect to the Weyl quantization: 

{Op]^ (a) f\g) = (nh)-'/' jj W{f ® g){V2x, -V20a{x, dx 

Even so, the same arguments give the same result, up to some minor re-scaling. Let 
be the operator given by 

M^{x,0 = {V2x,-V20. 

Then we have 

{U:W)if(E)g) = W{Utf®m^) 

= W{f®g) o o K^, o + 0{e)\\f\\L2\\g\\L^ 

= W{f^g)o~K.t + 0{e)\\f\\LA\9\\L^- 

Here 

kt := M^o Kt o M^, 
which we note is the Hamilton flow of the function 

P = -2poM^. 

Up to a re-scaling of t, we then have the Hamilton flow computed earlier, but now in the 
case = 8. 

We are most interested in HG modes and LG modes, so we take f ®'g = hmn- We then 
summarize our work in the following theorem: 
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Theorem 1. For T = -I'^l'^lt'l'^f = -2'^^'^h^^'^{a) + a - a)a)a - a)aa), having Weyl 
symbol p{x,$,; h) = \x{x'^ + — 2hx, and for Ut = e"**'^/'*, we have 

Here kt is the Hamilton flow of the function p given by 

p{x,^;h) = -2p{x/V2, -^/V2; h) 

+ - Ahx . 

Hence the intensities of the LG modes evolve along elliptic curves, as pictured in Fig- 
ure [21 This was in fact the motivation for the present paper. Previously, in [6j, we 
found that the Wigner transforms of HG modes evolve, in four-dimensional phase space 
(a;, ^, rf) e M^, according to a slightly more complicated flow, where the flow in the (a;, ^) 
plane (with x and ^ dual variables) for certain values of (y, Vj) appears as in Figure [2l Now 
we have found a situation where the flow along elliptic curves appears in the physical (x, y) 
plane. If one wishes, one may then consider Wigner transforms of LG modes; formulas 
for these are given in [5|. 
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